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METHOD AND APPARATUS TO GENERATE AN RF EXCITATION 
CONSISTENT WITH A DESIRED EXCITATION PROFILE USING A 

TRANSMIT COIL ARRAY 

BACKGROUND OF THE INVENTION 

The present invention relates generally to MR imaging and, more particularly, to a 
method and apparatus of parallel excitation by a transmit coil array to realize a desired 
excitation profile. The present invention further relates to a parallel excitation pulse 
design method that accounts for mutual coupling between coils of the coil array and 
applies to any coil geometry. The present invention is further directed to targeted RF 
excitation across an imaging volume to accelerate MR imaging. 

When a substance such as human tissue is subjected to a uniform magnetic field 
(polarizing field Bo ), the individual magnetic moments of the spins in the tissue attempt 
to align with this polarizing field, but precess about it in random order at their 
characteristic Larmor frequency. If the substance, or tissue, is subjected to a magnetic 
field (excitation field Bi) which is in the x-y plane and which is near the Larmor 
frequency, the net aligned moment, or "longitudinal magnetization", Mz, may be rotated, 
or "tipped", into the x-y plane to produce a net transverse magnetic moment M t . A signal 
is emitted by the excited spins after the excitation signal Bj is terminated and this signal 
may be received and processed to form an image. 

When utilizing these signals to produce images, magnetic field gradients (G x , G y , 
and G z ) are employed. Typically, the region to be imaged is scanned by a sequence of 
measurement cycles in which these gradients vary according to the particular localization 
method being used. The resulting set of received NMR signals are digitized and 
processed to reconstruct the image using one of many well known reconstruction 
techniques. 
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Spatially selective excitation is widely used in MR imaging to induce transverse 
magnetization while limiting the size of the signal-contributing volume. Slice-selective 
excitation, the most commonly used, confines the signal-contributing volume to a fixed 
slice that simplifies spatial encoding during signal acquisition to reduce data acquisition 
5 or scan time. Multi-dimensional excitation that produces localization along more than 
one dimension has been used to further this reduction in scan time. For example, 
localized spectroscopy, reduced-FOV scan of a region of interest, imaging of a target 
anatomy of unique shape, and echo planar imaging (EPI) with a shortened echo train 
length are applications usually implemented because of their support of reduced scan 
10 times. In addition, profile (flip, phase, and frequency) control across a sizeable volume 
with selective excitation has been exploited to improve excitation profile fidelity in the 
presence of Bo inhomogeneity or gradient non-linearity, and to reduce susceptibility 
artifacts. 



15 Selective excitation is commonly implemented with a single transmit coil that 

transmits across an entire volume and produces a relatively uniform Bi field, e.g., a 
birdcage coil. Highly efficient pulse algorithms have been developed for designing 
excitation pulses that suit such a configuration. Notwithstanding the advantages 
achieved by these pulse design tools, technical difficulties remain. Issues with excitation 

20 pulse duration, excitation profile accuracy, and RF power absorption (SAR) represent 
some of the outstanding challenges in a variety of applications. Compared to ID 
excitation, flexible profile control along multiple dimensions with 2D or 3D excitation 
entails intensified pulsing activity and often requires powerful gradients to keep pulse 
duration in check. This limitation hinders applications of multi-dimensional excitation on 

25 scanners with general-purpose gradients. Substantial subject-dependency of Bi field, 
resulting from increased wave behavior and source-subject interaction at high 
frequencies, may also contribute to the difficulty of excitation profile control. An 
elevated rate of RF power deposition at high frequencies represents yet another factor 
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that has a significant impact on the design and application of RF transmit modules and/or 
excitation pulses. 

Use of adiabatic pulses represents a pulse design approach that addresses the 
5 difficulty of excitation profile control associated with Bi inhomogeneity. This approach 
is limited as its application has been limited to certain profiles and tends to involve high 
RF power. A B\ -field optimization approach that aims at maximizing global B\ 
homogeneity addresses the control issue through transmit module improvements. 
Adaptation of the transmit coil geometry or the driving mechanism has been shown to 
10 reduce Bi inhomogeneity. At high frequencies however, the capability of a field 
optimization approach is limited. Even with calibration-guided adjustment of driving 
port weights, the degree to which the spatial variation of the composite Bi field 
approaches a desired level is highly dependent on the characteristics of component Bi 
fields, and results tend to be subject to considerable residual inhomogeneity. 

15 

Another proposed solution to reduce excitation pulse length is based on a parallel 
excitation architecture - multiple transmit elements driven by independent drivers. 
Individual Bi field patterns are employed to suppress aliasing lobes arising from 
sampling density reduction in the excitation k-space. Notwithstanding the excitation 

20 pulse length reduction achieved with a parallel excitation structure, application of such a 
structure has shown that particulars of the transmit elements are not fully taken into 
account. That is, these known parallel transmit architectures fail to consider mutual 
coupling between transmit elements and are often dependent upon a simplistic transmit 
array geometry. As such, spatial variations created by the transmit elements are not fully 

25 exploited. 

It would therefore be desirable to have a system and method capable of realizing 
desired excitation profiles and reducing excitation pulse length by the means of a parallel 
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transmit element architecture, where appropriate Bj field spatiotemporal variations are 
effected in a composite B ( field created by a transmit coil array. 
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BRIEF DESCRIPTION OF THE INVENTION 

The present invention provides a system and method of effecting spatiotemporal 
variations in a composite Bi field created by a transmit coil array that overcomes the 
aforementioned drawbacks. 

5 

The present invention is directed to the acceleration of a multi-dimensional 
excitation through the orchestrated driving of multiple transmit coils. The present 
invention emphasizes the coordination of multiple transmit elements to effect appropriate 
Bi spatiotemporal variations in a composite Bi field in order to improve the management 
10 of multi-dimensional pulse length while facilitating faithful production of desired 
excitation profiles. The present invention is also directed to the design of parallel 
excitation pulses with spatial and spatial-frequency domain weighting. 

Therefore, in accordance with one aspect of the present invention, a method is 
15 presented that includes the step of determining a desired RF excitation profile. The 
method further includes the step of independently driving each transmit coil of a transmit 
coil array such that the result of a collective excitation generated by the transmit coil 
array substantially matches the desired RF excitation profile. 

20 In accordance with another aspect of the invention, an MRI apparatus includes an 

MRI system having a magnet to impress a polarizing magnetic field and a plurality of 
gradient coils positioned about the bore of the magnet to impose a magnetic field 
gradient. An RF transceiver system and an RF switch are controlled by a pulse module 
to transmit and receive RF signals to and from an RF coil assembly to acquire MR 

25 images. The MRI apparatus also includes a computer programmed to design a plurality 
of RF pulse waveforms configured to control RF generation by a transmit coil array such 
that the result of collective RF generation across an imaging volume substantially 
matches a desired RF excitation profile independent of transmit coil array geometry. 
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In accordance with another aspect of the invention, the invention is embodied in a 
computer program stored on a computer readable storage medium and having instructions 
which, when executed by a computer, cause the computer to control RF transmission by a 
plurality of transmit coils of a transmit coil array such that spatial and temporal variation 
5 in a composite Bi field, accompanied by appropriate gradient changes played out in 
synchrony, create a desired excitation profile upon completion of RF transmission. 

Various other, features, objects and advantages of the present invention will be 
made apparent from the following detailed description and the drawings. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The drawings illustrate one preferred embodiment presently contemplated for 
carrying out the invention. 
In the drawings: 

5 Fig. 1 is a schematic block diagram of an MR imaging system for use with the 

present invention. 

Fig. 2 is a block diagram illustrating a linear transmit coil array assembly in 
accordance with one aspect of the present invention. 

Fig. 3 is a block diagram illustrating a wrap-around transmit coil array assembly 
10 in accordance with another aspect of the present invention. 

Fig. 4 is a graph illustrating an RF excitation profile achievable with a transmit 
coil array in accordance with the present invention. 

Figs. 5 - 6 are plots illustrating k x -direction weighting contribution by the coils of 
a transmit coil array positioned at two x-axis locations. 
15 Fig. 7 illustrates the magnitude of localization profiles along the x-axis for each 

coil of a transmit coil array. 

Fig. 8 graphically illustrates a pulse sequence in accordance with one aspect of the 
present invention. 

Fig. 9 illustrates resulting 2D transverse magnetization distribution as estimated 
20 by removing coil sensitivity weighting from an acquired image. 

Fig. 10 illustrates Bi field maps for the coils of an exemplary transmit coil array. 

Fig. 11 illustrates transverse magnetization distribution from a non-selective 
excitation in a reference body coil. 

Fig. 12 illustrates Bj field maps for each coil of a transmit coil array as well as a 
25 composite field map generated by superimposing the individual Bj field maps. 

Figs. 13-16 illustrate results of an RF pulsing protocol to control RF transmission 
and minimize RF deposition on a subject in accordance with another aspect of the present 
invention. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 

Referring to Fig. 1, the major components of a preferred magnetic resonance 
imaging (MRI) system 10 incorporating the present invention are shown. The operation 
of the system is controlled from an operator console 12 which includes a keyboard or 
5 other input device 13, a control panel 14, and a display screen 16. The console 12 
communicates through a link 18 with a separate computer system 20 that enables an 
operator to control the production and display of images on the display screen 16. The 
computer system 20 includes a number of modules which communicate with each other 
through a backplane 20a. These include an image processor module 22, a CPU module 

10 24 and a memory module 26, known in the art as a frame buffer for storing image data 
arrays. The computer system 20 is linked to disk storage 28 and tape drive 30 for storage 
of image data and programs, and communicates with a separate system control 32 
through a high speed serial link 34. The input device 13 can include a mouse, joystick, 
keyboard, track ball, touch activated screen, light wand, voice control, or any similar or 

1 5 equivalent input device, and may be used for interactive geometry prescription. 

The system control 32 includes a set of modules connected together by a 
backplane 32a. These include a CPU module 36 and a pulse generator module 38 which 
connects to the operator console 12 through a serial link 40. It is through link 40 that the 

20 system control 32 receives commands from the operator to indicate the scan sequence 
that is to be performed. The pulse generator module 38 operates the system components 
to carry out the desired scan sequence and produces data which indicates the timing, 
strength and shape of the RF pulses produced, and the timing and length of the data 
acquisition window. The pulse generator module 38 connects to a set of gradient 

25 amplifiers 42, to indicate the timing and shape of the gradient pulses that are produced 
during the scan. The pulse generator module 38 can also receive patient data from a 
physiological acquisition controller 44 that receives signals from a number of different 
sensors connected to the patient, such as ECG signals from electrodes attached to the 
patient. And finally, the pulse generator module 38 connects to a scan room interface 
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circuit 46 which receives signals from various sensors associated with the condition of 
the patient and the magnet system. It is also through the scan room interface circuit 46 
that a patient positioning system 48 receives commands to move the patient to the desired 
position for the scan. 

5 

The gradient waveforms produced by the pulse generator module 38 are applied 
to the gradient amplifier system 42 having G x , G y , and G z amplifiers. Each gradient 
amplifier excites a corresponding physical gradient coil in a gradient coil assembly 
generally designated 50 to produce the magnetic field gradients used for spatially 

10 encoding acquired signals. The gradient coil assembly 50 forms part of a magnet 
assembly 52 which includes a polarizing magnet 54 and a whole-body RF coil 56. A 
transceiver module 58 in the system control 32 produces pulses which are amplified by 
an RF amplifier 60 and coupled to the RF coil 56 by a transmit/receive switch 62. The 
resulting signals emitted by the excited nuclei in the patient may be sensed by the same 

15 RF coil 56 and coupled through the transmit/receive switch 62 to a preamplifier 64. The 
amplified MR signals are demodulated, filtered, and digitized in the receiver section of 
the transceiver 58. The transmit/receive switch 62 is controlled by a signal from the 
pulse generator module 38 to electrically connect the RF amplifier 60 to the coil 56 
during the transmit mode and to connect the preamplifier 64 to the coil 56 during the 

20 receive mode. The transmit/receive switch 62 can also enable a separate RF coil (for 
example, a surface coil) to be used in either the transmit or receive mode. 

The MR signals picked up by the RF coil 56 are digitized by the transceiver 
module 58 and transferred to a memory module 66 in the system control 32. A scan is 
25 complete when an array of raw k-space data has been acquired in the memory module 66. 
This raw k-space data is rearranged into separate k-space data arrays for each image to be 
reconstructed, and each of these is input to an array processor 68 which operates to 
Fourier transform the data into an array of image data. This image data is conveyed 
through the serial link 34 to the computer system 20 where it is stored in memory, such as 
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disk storage 28. In response to commands received from the operator console 12, this 
image data may be archived in long term storage, such as on the tape drive 30, or it may 
be further processed by the image processor 22 and conveyed to the operator console 12 
and presented on the display 16. 

5 

The present invention is directed to a method and system of accelerating RF pulse 
transmission by a plurality of transmit coils. Such a transmit coil array is illustrated in 
Fig. 2. Transmit coil array assembly 70 includes a plurality of RF coils or elements 72 
that are designed for parallel RF transmission, and a plurality of RF amplifiers 74. In one 

10 preferred embodiment, each transmit coil 72 is driven by a dedicated RF amplifier 74. In 
this regard, each RF amplifier is configured to generate a controlled current in a 
respective RF coil for defining and steering an excitation volume 76 of a subject 78 
within an MRI system. As will also be described, each of the transmit coils is controlled 
in a manner such that inter-coil correlations, i.e. mutual coupling, are taken into account. 

15 As illustrated in Fig. 2, the transmit coils 72 are arranged in a substantially linear fashion. 
Additionally, as will be described in greater detail, the RF amplifiers drive the plurality of 
RF transmit coils such that RF excitation may be substantially localized to a particular 
region-of-interest so as to reduce RF power deposition on the subject. As will be further 
described, each of the transmit coils is controlled in a manner such that RF power 

20 deposition is further reduced. 



Referring now to Fig. 3, transmit coil array assembly 70 is illustrated in another 
embodiment. In this embodiment, the transmit coils 72 are positioned in a wrap-around 
manner. In this regard, the coils are arranged in a distributed manner around the subject. 
25 Similar to that shown and described with respect to Fig. 2, each RF coil 72 is connected 
to a dedicated RF amplifier 74. One skilled in the art will readily appreciate that Figs 2-3 
illustrate a pair of possible arrangements of the coils of a transmit coil array and that 
other arrangements not specifically illustrated are possible and contemplated. 

10 
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As indicated above, the present invention is directed to a method and system 
operable with a transmit coil array such that RF excitation by the transmit coils is carried 
out in parallel. This parallel excitation supports a reduction in scan time through the 
5 acceleration of RF pulses and the localization of targeted excitation. 



The present invention will be described with respect to a small-tip-angle 
excitation, but one skilled in the art will appreciate that the present invention is 
extendable to other excitation regimes. The transverse magnetization resulting from a 
10 small-tip-angle excitation with a single transmit coil may be analyzed by the Fourier 
transform of the k-space trajectory traversed and weighted during the excitation: 

M(x) = jy M 0 (x) Z>(x)jV(k)S (k)e j2 " kx dk 

Eqn. 1, 

where S(k) represents a spatial-frequency sampling trajectory controlled by the switching 
15 gradients, W(k\ a spatial-frequency weighting induced by the driving RF source, and 
6(x), a spatial weighting induced by the coil's Bi field pattern. 



When several sets of pulse synthesizers and amplifiers form parallel RF sources 
that simultaneously drive corresponding coils during excitation, multiple spatial- 
20 frequency and spatial weightings influence the creation of the transverse magnetization. 
Within the limits of the small-tip-angle approximation, the k-space perspective expressed 
_ by Eqn. 1 may be extended to analyze a parallel excitation system based on the property 
of linearity: 

M(x) = jy M 0 (x) £ b n (x)\ ' f c fli/ ^(k)5(ky M ^ 

n=\ 1=1 

25 Eqn. 2. 
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In Eqn. 2, denotes the total number of transmit coils, n and / are coil indices, c n j are 
coefficients characterizing the mutual coupling between the coils, Wi(k) represent spatial- 
frequency weightings induced by the independently controlled RF sources, and b n (x) 
represent spatial weightings induced by the coils respective Bi field patterns. 

5 With g(x) denoting the term in Eqn. 2 that defines the excitation profile, g(x) may 

be expressed as 

*<*) = If wit ^(H^Ck^"^ = Z4(*)[ 

Eqn. 3, 

which indicates that in the analysis of the parallel transmit system, & / (x) = ^ ; c^A n (x)» 

10 the effective spatial weightings, may be used to account for coupling-induced inter-coil 
correlations. 



As an example, a 2D excitation case is considered, where an echo planar (k x ,k y ) 
trajectory, with k x being the slow direction and being the sampling period, is used and 
15 {(xy)\x min < x<x max , y m in<y<y m ax} specifies the field-of-view that contains the subject. The 
k-space weighting and sampling gives rise to a 2D excitation profile, which, as defined 
by Eqn. 3, is a weighted superposition of N periodic functions: 

g(x, y) = £ b l (x, y) £ u { (x - m A, y) Eqa 4 . 

/=/ m=-oo 

In Eqn. 4, the notation u^x) and A represent, respectively, J^/(k)e /27lk x dk and 1/A^. Z- 
20 dependence has been suppressed for simplicity. 

From Eqn. 4, it is clear that the discrete nature along k x necessarily implies 
aliasing lobes along x. Of significance, Eqn. 4 indicates that side lobe suppression may 
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be achieved through multiple weighting in the spatial (^(x)) and spatial-frequency 
(Wi(k)) domains. This can be compared to the case of excitation with a body-coil 
(volume coil with 6(x)«l), where a typical pulse design has the side lobes pushed outside 
the subject by limiting sampling period to be no greater than 1/D (D=x max -x min ). 

Within a small-tip-angle regime, design of gradient and RF pulses given a desired 
excitation profile may be achieved solving an inverse problem defined by Eqn. 3. For the 
purpose of illustration, a 2D excitation will be described. 



To achieve a 2D excitation profile given by g(xy) and with solutions of type: 
ufa>y) = hf(xy)g(xy\ Eqn. 4 may be rewritten as: 



+00 N 



g(*,y)= X^x-m^^^x-mA^foy) Eqn. 5, 

w=-oo l=\ 



which in general requires, for all (xy) inside the field-of-view, 

N 



Y J h l {x-m^y)b l {x,y)^ 



1, m-0 
0, otherwise 



Eqn. 6. 



By sorting the equations (e.g., through change of variables), it can be shown that 
{hfay), 1=1,. . .,7V} is typically constrained, at each (x^), by K linear equations (K is 
defined as the smallest integer that is greater or equal to D/A): 

C (x,y)\x,y) = e i Eqn. 7, 

where 
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b 2 (x,y) 



b N (x 9 y) 



b^x + mAty) b 2 (x + mA,y) ... ^(jc + mA,;;) 



Eqn. 8, 



h (x 9 y) =[ h \( x >y) h 2 (x,y) ... h N (x,y)] 7 

e l= [i o ... of 



Eqn. 9, 



Eqn. 10, 



and {x 9 x+mA (m*0), ...} represents the set of x coordinates within the field-of-view 
that are evenly spaced and inter-associated due to aliasing. Employing a sampling period 
&kx that is greater than 1/D, all but the first equation in Eqn. 7 represent the suppression 
of aliasing side lobes located within the field-of-view. 



Solving Eqn. 7 repeatedly for locations throughout the field-of-view yields 
10 /*/(xj>)'s, which then allow the calculation of k-space weighting according to the 
following: 

W, (k) = I h, (x)g(x)e^ k 'Vx Eqn. 11. 

The k-space weighting,, and the RF pulse waveform associated with the /th coil, can thus 
be calculated with the Fourier transform of a spatially-weighted version of the desired 
15 excitation profile, where the spatial weighting is derived from Bi field maps of each 
transmit coil and the k-space traversing trajectory. 



Quality of Bi field maps has a direct impact on excitation profile accuracy. The 
maps may be experimentally calibrated one at a time. With this approach, each 
20 calibration may involve an imaging experiment that uses a single element of the transmit 
array for transmission (with zero inputs to other elements) and the body coil for 
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reception. A division of the result by a reference image for removing the modulation of 
subject contrast and additional processing for suppressing the effects of noise, then 
provides an estimate of the effective Bi map associated with the transmit element. 
Alternatively, Bj maps may be inferred from sensitivity maps based on the principle of 
5 reciprocity. It should be noted that multiple sensitivity maps may be calibrated in 
parallel to reduce calibration time. However, the opposite phase and possible changes in 
coil coupling characteristics between transmit and receive, if not accounted for, may 
compromise the accuracy of the estimated effective Bi maps. 

10 Comparing two types of systems in the 2D excitation example, the present 

invention provides excitation acceleration of up to TV-fold over a single-channel body-coil 
system. Formally, this is revealed by the fact that Eqn. 7 admits at least one solution if 
N>D/A, or equivalently, A^^V/D, which is in contrast to the more stringent requirement 
of Aytc<l/D in the case of body-coil transmission. Intuitively, the capacity for 

15 acceleration, or, reduction in excitation k-space sampling density, is probably best 
appreciated by recognizing that while a reduction in excitation k-space sampling density 
causes aliasing lobes to locate inside the subject, an appropriate design of the spatial- 
frequency domain weighting (#7(k)) can combine with the spatial domain weighting 
(£,(x)) and the aliasing pattern (as determined by the sampling) to cause incoherent 

20 addition, therefore realizing reduction or annihilation of aliasing lobes' net amplitudes. 

For an acceleration factor that is smaller than N, or equivalently, a sampling 
period that is smaller than NfD, Eqn. 7 allows a family of solutions of dimensionality N- 
K. This results in choices of excitation pulse designs that are all capable of producing a 
25 main lobe that matches the desired excitation profile and, when applicable, 
simultaneously suppressing aliasing lobes. The specific design that uses A/(xj;)'s 
calculated by solving Eqn. 7 in the minimum norm sense is notable since it tends to 
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lessen the sensitivity of the excitation profile to perturbations or reduces the power 
requirement on the RF amplifiers. 

The independent driving of transmit coils of a transmit coil array also supports 
5 SAR management. Compared to uniform coverage of a subject volume with a single 
transmit coil, focused excitation of only the region-of-interest with an array of distributed 
local transmit coils by employing the coils in close proximity prevents substantial RF 
power deposition beyond the region. In addition, from the many ways of orchestrating 
the sources and achieving a desired excitation profile, the one that induces an E field with 
10 as small as possible an ensuing RF power deposition can be chosen. 

While the present invention supports a number of SAR reduction techniques, i.e. 
focused RF excitation, SAR management with a focus on the minimization of SAR 
averaged over the subject volume and the excitation period, which is defined by: 



15 will be hereinafter described in greater detail. In Eqn. 12, cr denotes tissue conductivity; 
p, density; V 9 the size of the irradiated subject volume; and P, the total number of time 
points used to quantify the temporal average. 



20 conducting material. At low frequencies, the fields inside the slab tend to be dominated 
by the incident fields, which are produced by the currents in the coils. Following a quasi- 
static approach in analyzing electric and magnetic near-fields, the fields may be 
characterized with a vector potential A: 




Eqn. 12, 



Given, for example, multiple loop coils placed facing the surface of a large slab of 
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Eqn. 13, 



10 



15 



where the line integrals over the currents in the coils are based on filament approximation 
of the coil conductors, and the fields are related to A by B=VxA and E=-dA/dt. In this 
case, the \E(xj?Af)\ 2 term in Eqn. 12 may be evaluated as: 



N 



H 



Eqn. 14, 



which is a quadratic form in [Ij(pAt) hip At) ... IrfpAt)], a vector with values of the 
current waveforms at time pAt. Sorting out the volume integral and temporal summation, 
SAR ave may be expressed as a quadratic function in the samples of the current 
waveforms: 



SAR =S*F S 



ave 



Eqn. 15, 



where superscript //denotes conjugate transpose, matrix F carries entries evaluated based 
on Eqns. 12 and 14, and vector s collects in a corresponding order a total of NxP samples 
of the current waveforms. 



Provided that the electric field scales linearly with applied source functions, a 
quadratic relationship in the form of Eqn. 5 between average SAR and source function 
samples generally holds. In the presence of biological objects or at high frequencies 
however, solving Maxwell's equations is difficult and construction of the F matrix may 
need to rely on calibration results or direct E field measurements. 



20 
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Given the dependencies of the absorption rate and transverse magnetization on the 
applied source functions, the determination of a set of coordinated source functions that 
produces the desired excitation profile while inducing minimum SAR is possible. In the 
small tip angle regime or its extension where a linear treatment of the Bloch equations is 
appropriate, closed- form solution exists for multi-dimensional excitation design, which 
obviates the task of searching a vast design space. 



Continuing with the previously described 2D excitation example, equations of the 
form of Eqn. 7, which stem from the requirement of creating the desired main lobe in the 
10 subject while avoiding aliasing lobes, collectively constrain the spatial patterns of A/(x)'s. 
Pooling these equations together thus gives the design constraints, which, in a matrix 
form, may be expressed as: 

C aiiKi = e *// E q n - !6. 

In Eqn. 16, C a ii is a block-diagonal matrix with C^/s on the diagonal and zeros 
15 everywhere else, and h a u and e a ii are vectors representing, respectively, concatenated 
h^'s and efs. If a moving sample of the weighting functions is carried out at a constant 
rate, the ?Tj(k(t))'s are proportional to the current waveforms. The Fourier transform 
relationship between the fT/(k)'s and the /*/(x)'s allows rewriting Eqn. 15 in terms of h fl //: 

SAR ave =h2/Vh a// Eqn. 17. 

20 The quadratic form remains as Fourier transform defines a linear mapping from hj{\) to 
W[(k). A variable sample rate would only modify entries of matrix V to match gradient 
amplitude changes. As such, pulse design for SAR management may be achieved by 
minimizing a quadratic function subject to a linear constraint: 

minimize h" Vh „ 

a " a " Eqn. 18, 

subject to C all h all = e all 

25 which may be solved using well-known numerical techniques. 
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Design principles for small-tip-angle parallel excitation pulses such as that 
described above were evaluated in simulation and phantom experiments. To evaluate the 
design principle for accelerated multi-dimensional excitation, parallel excitation with a 
5 transmit coil array was first examined in a simulation study. The transmit array was 
comprised of nine identical 19.8cm x 6.4cm loop coils that were placed on a flat form and 
lined up along the x-direction. This array faced a thin slab object below the array surface. 
2D excitation with a desired excitation profile across the object in the form of 
gi x )=gx(x)-g z (z) was approached with parallel excitation pulses. In this case, use of an 
10 echo planar k x -k z trajectory consisting of k x =constant lines evenly spaced by 4^, the 
negligible y- and z-direction B\ variation in the localized volume, and the separability of 
g(x) yielded solutions to Eqn. 1 1 of the form W^V)^Ukx,i{kx) -Ukz(k z ), where 

U^(k x ) = lh l (x)g x (x)e- J2 ^dx 
U kz {k z ) = \ z g z {z)e J2 ^dz 

For purposes of this first experiment equations of form Eqn. 7 were constructed and 
15 weightings over k x -k z were determined. RF pulse waveforms were then calculated based 
on Eqn. 11. As a reference, body-coil excitation pulses aimed at the same 2D localization 
were designed. 



The design principle for accelerated excitation was further evaluated in a phantom 
20 study, which was carried out on a 1.5 Tesla MRI scanner (CVi, GE Medical Systems, 
Waukesha, WI) with a setup very similar to that of the simulation study noted above. 
The transmit coil array of interest was of the same geometry and placed 3 cm above a 
water-filled 41x19x1 cm brick phantom. As the scanner only supported single-channel 
RF pulse transmission, the study examined parallel excitation indirectly, by mimicking 
25 simultaneous driving of the nine array elements through a series of nine single-channel 
experiments. Validity of the approach is ensured by the property of linearity in the small- 
tip-angle regime, which allows the prediction of the result of a parallel excitation 
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experiment from the superposition of transverse magnetization distributions observed 
from single-channel excitation experiments. 

Specifically, a single transmit/receive loop coil of size 19.8cm x 6.4cm was 
5 attached to the scanner's RF interface. During the nine experiments, the coil was placed 
and driven one configuration at a time, each with a position and RF pulse corresponding 
to one of the nine elements on the virtual coil array that were desired to simulate. After 
completion of every transmission, the coil was immediately switched to the receive 
function, whereas throughout the experiments the scanner's body coil was kept detuned. 

10 2D excitation and acquisition were carried out with a gradient echo sequence. From one 
experiment to another, excitation k-space traversing was kept the same (i.e., echo planar 
k x -k z trajectory with k x being the slow direction) but the weighting (RF pulse) was 
changed according to the excitation pulse design. 2D acquisition produced images that 
mapped out the water phantom along the x and z directions (and projected along y, the 

15 normal direction of the 1cm slab). 2D transverse magnetization distributions were 
quantified by removing the coils' sensitivity profiles from the images. The distributions 
were then superimposed to provide an estimate of the distribution resulting from the 
corresponding parallel excitation experiment. By the design of the study, coil coupling is 
not a factor. Bi maps that were estimated based on Biot-Savart Law were used in both 

20 the RF pulse calculations and the sensitivity profile removal. 

In another study on excitation acceleration, an all-around array geometry was 
examined. The array consisted of seven transmit elements that were distributed 
azimuthally on a wrap-around form inside a scanner's patient bore. Computer 
25 simulations evaluated 2D excitation designs that localize along both x and y dimensions. 
Coupling between elements was not negligible and was taken into account with a 
coupling matrix determined from mutual inductance calculations. The designs used the 
original Eqns. 7 and 1 1 . 
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Effectiveness of the SAR management scheme described previously as integrated 
in the parallel pulse design was further evaluated. The evaluation was carried out in the 
same fashion as the first simulation study except for the application of parallel excitation 
5 pulses of design type defined by Eqn. 18 instead of Eqn. 7. With the calculated A/(jc,z)'s, 
Eqn. 1 1 gave weightings over k x -k z , which in turn determined RF pulse waveforms. The 
resulting excitation profile and average SAR were compared to that of the first simulation 
study. 

10 A discussion of the results of the above-described experiments follow. Focused 

excitation of a 5cm by 5cm region centered at x=8cm and z^O inside the slab object was 
investigated in the first simulation study. Based on a body transmit coil, a reference- 
design employed pulses that traversed 57 k x =constant lines at 4^=1/31.6 cycles/cm. The 
x-direction localization that resulted from this reference design is shown in Figs. 4-7. A 

15 parallel excitation design accomplished the 2D localization task with the transmit coil 
array. Representing a 4-fold acceleration, the design employed pulses that traversed 14 
k x =constant lines at 4^=1/7 cycles/cm. UkxAmAkx) and Ukxj(rnAkx\ the /^-direction 
weighting contributed by the coils positioned at x=-4cm and x=8cm, respectively, are 
illustrated in Fig. 5 and Fig. 6. Localization along x due to each of the nine coils is 

20 shown in Fig. 7. Note that while the first aliasing side lobes were 4.5 times closer to the 
target (center-to-center spacing = 7cm) as a result of the sampling density reduction, the 
net amplitudes of these as well as other aliasing lobes located inside the 40cm FOV were 
negligible due to incoherent addition, as shown in Fig. 4. Compared to the result of the 
body-coil approach, localization of the parallel excitation was as well refocused (the 

25 imaginary component, not shown, was negligible) and of comparable spatial resolution. 
See Fig. 4. 
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In the phantom study, effects of incoherent addition on aliasing side lobes were 
the focus of investigation. To this end, 2D excitation pulses were designed to target a 
region in the water phantom directly below the center element. To facilitate the 
investigation, pulse calculations further assumed an extended linear array instead of the 
5 9-element one. The designed pulses were 5.7msec in length. For the center element 
experiment, Fig. 8 shows the applied RF pulse (magnitude and phase) as well as G x and 
G z , the gradient pulses identically executed in all the experiments of the series. 
Removing the coil's sensitivity profile from the resulting image provided an estimate of 
the 2D transverse magnetization distribution induced by the element, as shown in Fig. 9. 

10 Fig. 10 illustrates the Bi /sensitivity maps used. As a reference, Fig. 11 illustrates the 
transverse magnetization distribution from a nonselective excitation in a body-coil 
transmit-receive experiment. Noticeable in Fig. 9 is a noise amplification effect due to 
the division operation employed for sensitivity profile removal, which tends to increase 
in severity farther away from the sensitive region. To prevent excessive noise 

15 amplification from obscuring the investigation, the division operation was suppressed in 
distant regions. 

Results from all nine experiments are summarized in Fig. 12, which displays in 
rows 1 through 9 the mapped transverse magnetization corresponding to each of the 

20 experiments. The bottom row (row 10) presents the result of superimposing the 
individual maps, intended as a prediction of the result of a corresponding parallel 
excitation. Again, substantial reduction of aliasing side lobes due to incoherent addition 
was observed. With the setup, contributions from the elements in the establishment of 
the main lobe and the suppression of the aliasing lobes were readily appreciated. The 

25 results from the center element alone and from the middle five and middle nine elements, 
suggest that local excitation profile control is mainly achieved through nearby coils. Use 
of the extended array assumption in the pulse calculations accounted for much of the 
residual , aliasing (incomplete annihilation) towards the 9-element array's boundary. 
Augmenting the array with elements beyond the nine can rectify this effect. Designing 
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pulses for the 9-element array can rid this effect too, in which case boundary coils' 
weighting would experience the greatest changes. 

2D parallel excitation pulses for a wrap-around array were designed and 
5 evaluated. The simulations concentrated on the task of selectively exciting an arbitrarily 
positioned local volume within a 40cm-by-23cm axial field-of-view. Eqn. 7 was solved 
repeatedly based on the effective B\ field patterns and an EPI trajectory comprising 14 
k x =constant lines at 4^=1/6.9 cycles/cm. For the /th coil, /=1,2,.. .,7, the product of the 
desired 2D localization profile with the calculated hj(xy) was then Fourier transformed to 
10 derive the coil's k-space weighting and RF pulse waveform by the parallel excitation. 
The net result was substantially free of aliasing side lobes and represents an excellent 
match to that of a reference excitation, which involved body-coil transmission of a 4- 
times longer conventional RF pulse. 

15 The design of the last simulation study resulted in parallel excitation pulses that 

differed in shape from the pulses of the first simulation study. Figs. 13-16 present the 
outcome with a format similar to that of Figs. 4-7. While the pulses maintained the same 
level of localization accuracy and spatial resolution as that of the pulses of the first 
simulation study, Fig. 13, the design changes led to a 38% reduction in average SAR, 

20 confirming the substantial impact of the integrated SAR management scheme. 

With the present invention, designed RF pulses are synthesized, amplified and fed 
to corresponding transmit elements in parallel to induce both spatial and temporal 
variations of the composite B\ field, which, accompanied by appropriate gradient changes 
25 played out in synchrony, create a desired excitation profile upon completion of excitation. 
This is in contrast to a conventional approach, where the design of coil geometry and the 
offsets of driving-port phase/magnitude target Bi -field spatial homogeneity, and an RF 
pulse played during excitation is limited to manipulate Bj -field temporal variation only. 
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One skilled in the art will recognize that inducing appropriate Bi spatiotemporal 
variations for excitation bears significant ramifications on RF excitation performance. 
That is, parallel excitation accommodates excitation acceleration and/or SAR control 
without substantial sacrifice in the accuracy of producing the desired excitation profile. 

5 

In summary, the RF pulse driving a transmit element can be calculated with the 
Fourier transform of a spatially weighted version of the desired excitation profile, the 
capacity for accelerating multi-dimensional excitation by the means of k- space sampling 
density reduction lies with the suppression of aliasing lobes and can be achieved by 
10 appropriately designed driving pulses (spatial-frequency domain weightings), and SAR 
management can be accomplished by minimizing a quadratic function in the driving 
sources, which searches a way of orchestrating the sources to achieve a desired excitation 
profile and/or acceleration while inducing an E field with minimum ensuing RF power 
deposition. 

15 

From an application perspective, fast imaging is an area where the present parallel 
excitation approach is particularly applicable. Under circumstances where the anatomy 
of interest is contained in a local region for example, multi-dimensional excitation that 
"spotlights" the region allows acceleration of imaging by alleviating the burden of spatial 

20 encoding inflicted on signal acquisition. Representing improvements over conventional 
excitations, multi-fold shorter parallel excitations support imaging volume 
definition/steering while breaking the time cost barrier that hindered the practical use of 
multi-dimensional pulses in the past. Compared to the use of a parallel acquisition 
approach, focused imaging based on the parallel excitation approach is not subject to the 

25 unique SNR degradation described by the geometric factor. Combined use of the two 
approaches is possible and can provide an even greater capacity for scan time reduction. 
While the experiments reported here focused on 2D localization, the parallel excitation 
approach applies to the creation and acceleration of general 2D excitation profiles, with 
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utilities including correction for field imperfection-induced effects and non-Fourier 
spatial encoding. The present invention is also applicable to 3D excitations. 

In high field imaging, the transmit system and driving means described may also 
5 be used to both manage excitation profile and regulate RF power deposition. Embodying 
an integrated treatment of excitation pulses and transmit coils, the present invention 
facilitates excitation profile control. Transmission with a distributed parallel system, 
acceleration of excitation and management of SAR further provides a solution to power 
deposition at high field strength. 

10 

The present invention further provides a method of MR imaging that includes 
determining a desired RF excitation profile and independently driving each transmit coil 
of a transmit coil array such that the result of a collective excitation generated by the 
transmit coil array substantially matches the desired RF excitation profile. As described 

15 with respect to Figs. 2 and 3, each transmit coil is connected to a dedicated RF amplifier 
that is designed to control RF excitation of the corresponding transmit coil. More 
particularly, each RF amplifier or other control means provides a control signal 
representing an RF pulse waveform specific to a corresponding transmit coil such that a 
composite RF excitation is generated across an imaging volume that coincides with the 

20 desired RF excitation profile. The method further includes determining the RF pulse 
waveform for each transmit coil from a spatially weighted version of the desired RF 
excitation profile. As mentioned previously, the spatial weighting of each transmit coil is 
considered such that mutual coupling between the coils of the transmit coil array is taken 
into account. 

25 

The present invention is also directed to an MRI apparatus that includes an MRI 
system having a magnet to impress a polarizing magnetic field, a plurality of gradient 
coils positioned about the bore of the magnet to impose a magnetic field gradient, and a 
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RF transceiver system and an RF switch controlled by a pulse module to transmit RF 
signals to an RF coil assembly to acquire MR images. The MRI apparatus also includes a 
computer programmed to design a plurality of RF pulse waveforms configured to control 
RF generation by a transmit coil array such that the result of collective RF generation 
5 across an imaging volume substantially matches a desired RF excitation profile 
independent of transmit coil geometry. The computer is also programmed to determine 
an effective spatial weighting imposed on the collective RF generation by each transmit 
coil of the transmit coil array and apply the plurality of RF pulse waveforms such that 
coupling induced inter-coil correlations are taken into account. The effective spatial 

10 weighting imposed by each transmit coil includes spatial-frequency weightings induced 
by RF sources configured to control RF excitation of the transmit coils, weighting 
associated with mutual coupling between the transmit coils of the transmit coil array, as 
well as a respective Bi field for each transmit coil of the array. The computer is also 
programmed to design the plurality of RF pulse waveforms given any transmit coil array, 

1 5 as such, the transmit coils may be linearly arranged, arranged in a wrap-around fashion, 
or any other arrangement that may be implemented. Additionally, one or more of the 
transmit coils may also be designed to receive MR signals. 



The present invention may also be embodied in a computer readable storage 
20 medium having a computer program stored thereon and representing a set of instructions 
that when executed by a computer causes the computer to control RF transmission by a 
plurality of transmit coils of a transmit coil array such that spatial and temporal variation 
in a composite Bj field induces a desired excitation profile upon completion of RF 
transmission. The set of instructions further causes the computer to control application of 
25 gradients in an imaging volume to be in synchrony with spatial and temporal variation 
creation in the RF transmissions. The set of instructions further causes the computer to 
control application of control signals to the plurality of transmit coils such that RF 
excitation by the plurality of transmit coils occurs in parallel. The computer is also 
caused to determine the spatial and temporal variations to be induced with at least one 
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transmit coil's effective Bi field map. The effective Bj field maps for each transmit coil 
reflect mutual coupling of the plurality of coils. Preferably, the effective Bj field maps 
are generated during calibration of the transmit coil array. 



5 The present invention has been described in terms of the preferred embodiment, 

and it is recognized that equivalents, alternatives, and modifications, aside from those 
expressly stated, are possible and within the scope of the appending claims. 
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